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The phase diagram of doped manganate compounds Lai-^A^MnOs (with divalent A) is studied. 
We analyze an extension of the double exchange model using the Schwinger boson formalism. Earlier 
work by de Gennes on the existence of a canted phase is reproduced, although this phase is shown to 
be unstable towards phase separation in a broad regime of physical interest. We numerically solve 
the mean field equations for our model and exhibit its phase diagrams. 
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Doped manganese oxides show many unusual features, 
the most striking being the colossal magnetoresistance 
in the ferromagnetic phase |l] || . The phase diagram, 
as function of doping and temperature is far from elu- 
cidated. At small dopings, many experiments are inter- 
preted in terms of the phase diagram proposed by de 
Gennes jQ, who studied the so called double exchange 
model |5[| (see below). Some experiments indeed confirm 
the predictions derived from this approach ||. Others, 
however, seem to imply a more complex behavior, includ- 
ing charge ordering ]7],|j or coexisting phases @JlQ]. In 
addition, these materials show a metal-insulator transi- 
tion at low dopings and low temperatures . 

In the following, we will analyze the phase diagram of 
these systems using the Schwinger boson representation 
p2[ for the magnetic moments, which are described by 
the double exchange model. We neglect the role of lat- 
tice distortions, which may be important at high dopings, 
where the Jahn- Teller distortion present in undoped sys- 
tems disappears JDjfl . The scheme that we use allows us 
to obtain a description of the spin waves. In more con- 
ventional systems, it has been shown that quantum and 
thermal fluctuations are adequately described |l2|]l4| l in 
this approach. The method has already been used to 
study the quasiparticle coherence in the manganese ox- 
ides jlB). Finally, the calculations reported here are in 
general agreement with the work of de Gennes [Q in the 
in zero temperature, large S limit. The general features 
of the model are described in the next section. Sections 
III and IV adapt the Schwinger boson method to the 
double exchange model. The results for T = are pre- 
sented in section V, where the relation of our work to the 
original calculation by de Gennes Q is discussed. Sec- 
tion VI is devoted to finite temperature results. Finally, 
section VII contains a discussion of experimental results 
and related thorctical work. 



In materials such as Lai-^A-rMnOs (with A divalent), 
a fraction x of Mn ions are in 3d 4 (Mn 3+ ) configura- 
tions, with the remaining fraction 1 — x in 3d 3 (Mn 4+ ) 
states. In a cubic crystal field, the Mn 3d levels split into 
a lower ti g triplet and an upper e g doublet. Intra-atomic 
("Hund's rules") couplings overwhelm the crystal field 
splitting, hence the ti g levels are always triply occupied 
and form a S — | 'core spin'. In Mn 3+ ions, the e g or- 
bitals are further split by a static Jahn- Teller (JT) distor- 
tion, which, together with the Hund's rules, completely 
determines both the orbital as well as spin state of the e g 
electron. The e g electrons may be represented by spin- 
less, single-orbital fermions whose hopping is modulated 
by the overlap of the core spin wavefunctions. If we treat 
the core spin on site i using the Schwinger representation, 
S t = ¥l° a php (M=T4, EcAU Q = 25), then the 
e g . electron creation operator if>j may be factored into 
a spinless fermion c and the Schwinger boson b ia which 
supplies the core spin orientation: ipj a = c\b iry . The role 
of the core spin overlap to electron hopping in these ma- 
terials is widely appreciated (see e.g. refs. [p|, ^3|JT^ , p^ ] ) . 

Neighboring core spins are coupled via superexchange 
through the O 2p orbitals @|T§. For pure LaMn0 3 
(x = 0), the c-axis exchange is antiferromagnetic whilst 
the exchange between neighboring ions in a plane perpen- 
dicular to c is ferromagnetic. We have therefore chosen 
to study the model defined by the Hamiltonian |18[ 



H.c 



Si • Sj 



for a cubic lattice of Mn ions, where the exchange Jy = 
— J v < along vertical links and Jy = Jh > along 
horizontal links. Note that this model is somewhat unre- 
alistic in that the exchange between core spins is fixed in- 
dependent of the fermion occupancy and hence it cannot 
reflect the difference between Mn 3+ -Mn 3+ , Mn 3+ -Mn 4+ , 
and Mn 4+ -Mn 4+ exchange (note that Mn 4+ -Mn 4+ ex- 
change, appropriate to pure CaMn03, is always antifer- 
romagnetic In addition, we assume a strong static 
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JT distortion which renders the conduction orbital non- 
degenerate, whereas in the real materials this distortion 
vanishes for x > 0.2. Nonetheless, the model does capture 
what is perhaps the most important aspect of the inter- 
action between fermions and core spins, namely that fer- 
romagnetic core spin alignment leads to a larger fermion 
bandwidth and reduced kinetic energy. We aim to ap- 
ply it in the small x region of the phase diagram, where 
most of the links are between Mn 3+ ions. This hopping 
Hamiltonian itself, in the absence of Heisenberg exchange 
terms, was considered by Sarker JL5| , who found a finite 
temperature transition between a ferromagnetic metal 
and a spin-disordered state, presumably insulating, in 
which the fermion band is completely incoherent. 

III. MEAN FIELD THEORY 

Following Sarker , we invoke a Hartree-Fock decou- 
pling of the hopping term: 

4 c j^ij -(4 c j)(Fii) + c l c j( :F v) + (cjc^Tij +flucts. , 

where T L j = &, tr ^i cr accounts for the overlap of the core 
spin wave functions. The Heisenberg exchange is treated 
using the static Schwinger boson mean field theory of ref. 
p2| . The Hartree-Fock Hamiltonian is then given by 

Tt HF = Ho + Whop + Wbosc + T~i cond 

with 

H = NS 2 J v + 2NS{S+l)J h -2S^2A t 

i 

+TE h i^fe-)i 2 + TE v iQvfo-)i 2 

+i5^[t ii (c]c.)(^)+H.c] 

m 

Whop = -\SY, c h (Fij) + H.c] 

(ij) 

Wbosc = ^[Qh{ij)^ij + H.c] + ^[Qvfa'Mij + H.c] 

(ij) {ij) 

+ E A > h lK - {c\ Cj ) + H.c] 

i,a (ij) 

k.cr 

Here, (ij) anc ^ (ij) represent sums over horizontal 
and vertical links, respectively. The are Lagrange mul- 
tipliers which enforce the local constraints b\ a b ia = 2S, 
Aij = {b^b-^ — bj^b i ^) measures the antiferromagnetic 
correlation between sites i and j, N is the total num- 
ber of sites, and B^a is a field which is conjugate to the 
Schwinger boson condensate order parameter: 



where F is the free energy. We are guided to a simple 
mean field theory with seven parameters, in which we 
assume 

^h,v = I*h,v^ ( C l c j)h,v 

Qh = Q(ij) on horizontal links 
Q v = e l7r Ri Q(ij) on vertical links 

are all real constants, where 7r = (0, 0, it) in units where 
the lattice constant is unity. Thus, there are seven mean 
field parameters. 

Diagonalizing 7Ybosci we find 

n Bosc = £ E(k) pij ku + J2 ( v A l - A l- A 

fc.cr k 



with 




and 

Afc = A — 2(Qh + ATh)(cos k x + cos k y ) 
Afc = — 2Q v cosfc z 
= —2K V cos k z 

E(k) ee y/A% -Al + n k . 

When there is a condensate (T < T c ) the Bose spectrum 
is gapless, with A = A*, where 

A* = 4(Q h + K h ) + 2^Ql + Kl . 

The dispersion then may be compared in the x — > limit 
with the spin wave result 

E sw (k) = S^[J V + J h (2 - cos k x - cos k y )} 2 - J 2 cos 2 k z ; 

obtained by expanding about a (0, 0, 7r) Neel state (alter- 
nating ferromagnetic planes). The basic functional de- 
pendence on k is reproduced - this is a good preliminary 
check on the mean field Ansatz. 

Note also the particle hole symmetry present in our 
mean field theory. This guarantees an x — > 1 — x symme- 
try in the phase diagram. As mentioned above, exchange 
in pure CaMnC>3 is different than in pure LaMnC>3, due to 
the presence of the second set of empty e g states. Hence, 
this symmetry is an artifact of our model. 
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In deriving the mean field equations for our model, 
we must also include the condensate. The relationship 
between the held B^ and the order parameter is 



7T-fcT 



M~\k) 



B 



7r-fc| 



B 



ki 



(2) 



Differentiating the condensate contribution to the free 
energy with respect to a generic mean held parameter £ 
gives 



dM(k) 



U„,b.,, it is 

(7(7 1(7 ' 



Enacting a global SU(2) rotation b ia . 
easy to show that the free energy is invariant under such 
a transformation. This approach to Schwinger boson con- 
densation can also be applied to the cases of the uniform 
ferro- or antiferromagnet. It has the comforting feature 
of making the SU(2)-invariance manifest from the out- 
set (compare e.g. with ref. |l4j|, in which the condensate 
always results in a moment in the x-direction). 

Proceeding in our analysis, we assume condensation 
only at k = and k = n. In order that the condensate 
give no contribution to the free energy, we require that 



7Tl 



01 



-le' 7 



-Ye 1 



cos 

- sin \d 

cos i?9 

— sin 



(3) 



where tan?9 = Q v /K v and X, Y, 7, and 6 are at this 
point arbitrary parameters specifying the direction and 
magnitude of what is in general a canted (0, 0, n) antifer- 
romagnet Q . Equation p| also is consistent with the free 
energy being a convex function of the order parameter 
^ka- The condensate is then spatially varying, with 



(4) 



The condensate contribution to the local magnetization 
(Si) is easily computed and $ is found to be the canting 
angle. 



IV. MEAN FIELD EQUATIONS 

We now are in a position to write down the mean held 
equations. We work in the grand canonical ensemble, 
introducing a chemical potential fi for the fermions. This 
introduces an eighth parameter. However, we hnd that 
the mean held equations guarantee Qh = ^Jh^h always, 
so we are left with the following seven equations: 
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d 3 k 



A, 



(27T) 



ctnhfffi) 
\2k n TJ 



R z 



T v = 2 



d 3 k Q v cos 2 k z 



ctnh 



(2tt) 3 

d 3 k (cos k x + cos k y ) A fc 
(2^)3 

d 3 k 



cos k. 



( E(k) \ 
\2k B TJ 

ctnh ( 



Q V R 2 



E(k)- 



2k B T 
K V R 2 



R' 



(2tt) 3 



Kh = jSt h 



cx p(I^) 

d 3 k cos k r 



1 VQ$ 

cos ky 



(27T) 



\ k B T J 



exp( 

d 3 k cos k 7 
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(27T) 



d 3 k 



exp 
I 



\ k B T J 
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exp 
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where R = \/X 2 + Y 2 is the condensate amplitude, x is 
the hole concentration, and 

e(fe) = — S , ^h-? r h(cos k x + cos k y ) — St v T v cos k z 

is the fermion dispersion. The integrals are performed 
over the hrst Brillouin zone of the cubic lattice. There 
are seven mean held equations corresponding to seven 
mean held parameters. The parameters are: 

T<T C : Q v ,T h ,T v ,K h ,K Vlf i,R (A = A*) 
T>T C : Q v ,T h ,T Vl K h ,K v ,n,A (R = 0) . 

We have identihed several phases which emerge from 
the mean held theory: 

(I) Antiferromagnet (LRO at n) : T v , K v = 0, R ^ 
(II) Canted (LRO at and tt) : Q v , K V ,R^Q 

(III) Ferromagnet (LRO at 0) : Q v = 0, R ^ 

(IV) 3d Local Magnetic Order : F v , F h ^ 0, R = 
(V) 2d Local Magnetic Order :F V ,R = 0, -F h ^ 

(VI) Maximally Disordered : Q v , !F V , JFh, K v , K^, R = . 

The ordered phases I, II, and III were identihed by de- 
Gennes M; we have also found evidence of phase sep- 
aration below T c (sec also 19 20[]). In what follows we 
describe our analytical and numerical investigations of 
the phase diagram. 



V. T = 0, 5" -» 00 LIMIT 

Our mean held equations simplify considerably in the 
limit of zero temperature and S — > 00. We examine the 
three ordered (A = A*, R > 0) phases, 

(I) Q v = SJ V , T h = 2S, F v = K v = 
(II) Ql + Kl = S 2 J 2 , = 2K V /J V , JT h = 25 
(III) Q v = 0, T y = T h = 2S . 
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The canted phase II can smoothly interpolate between 
the 7r-antiferromagnet I and the ferromagnet III, with $ 
going from \-k to 0. We start with the canted structure, 
solving the mean field equation T v — 2K V /J V . We do this 
in the regime x *C 1 by expanding the fermion dispersion 
relation 

e(fe) =e(0) + St v T v {\- cos k z ) + ±St Yl T^{k 2 x + k 2 y ) + . . . : 

since T w = is a possible solution, we keep the full c-axis 
dispersion. We find that the solution is characterized by 
the dimensionless parameter r = 87rJ v ihA v - For r > 1, 
the only solution has T v — K v — 0, and we have antifer- 
romagnetism at finite doping. However, experiments ]22] ] 
suggest J v w 0.58 meV while spin density functional cal- 
culations j23| suggest t v w 44 meV (the physical hopping 
parameter is S 2 t w 100 meV). This gives r w 0.33, so 
the r > 1 regime is unphysical for Lai_ x A K Mn03. The 
ground state energy per site is found to be 

— ^ = -2J h - J v - 4£ h .x + 2nt h x 2 + ... (6) 

For r < 1, we have a solution with nonzero J- v . We 
find 

^ = -2 J h - J v - 4i u x + 2ng{r)t h x 2 + ... (7) 

7r 2 r — Tree sin 2 a 

ff( r ) = T- vT 

(sin a — a cos a) z 

where a(r) € [0,7r] is defined implicitly by the equation 

a — sin a cos a = nr . (8) 

We obtain T v = xt v /J v + 0{x 2 ) in this regime. For 
r < 5) <?( r ) < and the coefficient of a; 2 is negative, 
indicating phase separation. For i < r < 1, g(r) > 
and there is a homogeneous, thermodynamically stable 
canted phase, originally identified by deGennes Our 
estimate r w 0.33 suggests that phase separation is likely. 

Within the canted phase, as x increases from to \, the 
canting angle decreases from d — ^ir (7T-LRO) to d = 
(0-LRO). The transition from canted to ferromagnetic 
order is continuous, occuring at a critical concentration 
x* determined by the simultaneous solution of the two 
equations 

/rl 3 k 
jSjseic-'W) 

in the two variables x* and fi. For sufficiently large J v , 
is nonzero even at half-filling and the system remains 
canted for all x. 
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FIG. 1. Energy per site E/N (solid) versus concentration 
x for J v = Jh = J, t v = th = t = 1 at temperature T = 0.01 1 
for three different values of r = 8nJ/t: (a) r — 0.25, (b) 
r = 0.50, (c) r = 0.75. To ascertain the sign of d 2 E/dx 2 , we 
have subtracted from E/N the linear part; the dashed curve 
is the remaining contribution. Note that E(x) is convex for 
small x in (a) and (b), indicative of phase separation. 

VI. FINITE T PHASES 

To explore the finite temperature phases of our model, 
we have solved the mean field equations numerically us- 
ing the MINPACK routine hybrd.f. (To simplify matters, 
we assumed J v = Jh = J and t v = = t.) We found 
that there is often more than one solution to the mean 
field equations; in such cases we computed the energy of 
each solution and identified the minimum energy state. 
To compare with the analytical results of the previous 
section, we computed the energy versus concentration x 
for three different values of r, all at a temperature well 
below T c (see figure 0). This verified our prediction of 
phase separation for < r < |. 

In figures 2 and 3 we plot phase diagrams for J/t = 
0.01 (r = 0.251) and J/t = 0.03 (r = 0.754), respec- 
tively. Regions are labeled I through VI corresponding 
to the six phases discussed above (recall that a conden- 
sate is present only in phases I, II, and III). We find that 
only region III of figure 2 is stable with respect to phase 
separation. The homogeneous mean field solution yields 
a convex F(x, T) outside of this region up to tempera- 
tures on the order of T*{x) < 0.86 1; for T > T* the free 
energy is concave in x. 

We found a first order line separating phases I and II 
from the disordered phases. This is so even at x = 
the mean field theory predicts a first order transition from 
the (0, 0, 7r) Neel state to a magnetically disordered state. 
This is perhaps a worrisome artifact of the mean field the- 
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FIG. 2. Phase diagram for J/t = 0.01 (r = 0.251) ob- 
tained from numerical solution of the mean field equations. 
The dark solid line separating phases II and IV is first order. 
All other transitions are second order. Dotted lines repre- 
sent transitions between disordered states. Phase separation 
occurs outside region III (see text for discussion). 
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FIG. 3. Phase diagram for J/t = 0.03 (r = 0.754) ob- 
tained from numerical solution of the mean field equations. 
A small sliver of the 7T-ordered phase exists in the lower left 
corner. The dark solid line separating phases I and IV as well 
as II and IV is first order. All other transitions are second or- 
der. As in figure 2, dotted lines represent transitions between 
disordered states. 



ory. In addition, the transitions between the disordered 
phases may well become smooth crossovers when fluctu- 
ation effects are accounted for. Indeed, application of the 



Schwinger boson formalism to the Heisenberg model |lj] 
leads to a spurious high temperature mean field transi- 
tion to a state in which the magnon bandwidth vanishes, 
analogous to phase VI. 

Sarker has discussed the behavior of the electron 
spectral function and found it to be entirely incoherent 
when the core spins are disordered. Writing (at finite 
temperature T) 

G^(k,t) = -i({^ ka (t),^ kf3 (o)})e(t) 



du>' 



i0+ 



we find, at the mean field level, 
d 3 q 



(2tt) 3 



(q)+S a(3 f(k + q) 



xd (w + E{q) - e{k + q) + fi) , (10) 



where n a/3 (q) = (b ka b k/3 ) and f(p) = (cj,c p ) are equi- 
librium averages. The contribution of the condensate to 
the spectral density results in well-defined quasiparticle 
peaks. For instance, 

p™ nd {k, w) = Y 2 sin 2 (tf/2) S(u + E{0) - e(k) + (i) 

+X 2 cos 2 (z9/2) S(w + E(tv) - e(7r - k) + (J,) 

where X, Y, and $ describe the amplitude and orienta- 
tion of the condensate (recall equation]^). The remaining 
contribution to the spectral function, Ap, is incoherent 
and spectrally broad |l5j . Our calculation allows for con- 
densation both at k = as well as k — 7r, and as expected 
there are two quasiparticle peaks when translational sym- 
metry is broken (phase I). A detailed study of p(k,ui) in 
the various ordered and disordered phases is pending. 



VII. CONCLUSIONS 

We have shown that the double exchange model has 
a variety of possible phase diagrams, controlled by the 
parameter r = QirJ^th/t 2 .. For realistic values of r, 
0.05 < r < 0.2, we find a phase diagram similar to the one 
proposed by de Gennes , except that the canted phase 
is replaced by a region of phase separation. In addition, 
the transition to the paramagnetic phase may be of first 
order, and the paramagnetic phase itself is anisotropic. 
For r > i, we find that the canted phase is stable. The 
main source of uncertainty in the value of r arises from 
the lack of a precise determination of the hoppings, which 
may depend on the composition and details of the lattice 
structure j3|. Perhaps both situations may be realized 
experimentally. 

Coulomb interactions will prevent charge separation at 
large scales. The electrostatic energy required to break 
the system into charged domains of side t goes as e 2 x 2 /e£, 
where e is the dielectric constant. The magnetic energy 
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cost to create a domain wall of size t is roughly J v (£/a) 2 , 
where a is the lattice spacing. Hence, the domain size, £ 
goes as a 4 / 3 (J v e/e 2 x 2 ) 1 / 3 , which should be on the order 
of a few lattice constants. 

Phase separation in these compounds has previous 
been discussed phenomenologically in jl9), and in |?o|| , 
in the context of numerical results for the ferromagnetic 
Kondo lattice in one, two, and infinite dimensions. At 
large values of the Hund's rule coupling, this model re- 
duces to the double exchange model, plus antiferromag- 
netic interactions between the core spins. In-plane ferro- 
magnetic interactions do not arise, as they are induced by 
the second e g band. Our results, for T — 0, are qualita- 
tively in agreement with those reported in |^0| , provided 
that one identifies our canted phase with the incommen- 
surate order reported there. 

Our calculation reproduces the observed magnon spec- 
trum at zero doping pi] , except for a small anisotropy 
gap. In the canted phase, the long wavelength spin waves 

behave as J {k 2 + k 2 ) + v ± k 2 . In the phase separated 

regime, localized ferro- and antiferromagnetic modes are 
expected, which may have been observed experimentally 

A detailed study of transport properties lies beyond 
the scope of the present calculation. It is interesting to 
note, however, that, in the presence of phase separation, 
hopping is suppressed in the out of plane direction in the 
antiferromagnetic domains. This effect can contribute 
to make these compounds insulating at low dopings, in 
agreement with experiments. 
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